function a = coef(x, y)
	% Intrari:
	% x = tabloul absciselor celor n puncte
	% y = tabloul ordonatelor celor n puncte
	% Iesiri:
	% a = coeficienti polinom Lagrange
	n = length( x );
	a = zeros( n, 1 );
 	z = zeros( n, 1 );
	 % calcul coeficienţi c din (x-x(1))...(x-x(n))
	c = poly( x );
 	for i = 1 : n
		% calcul coeficienţi b ai împărţirii
 		% polinomului prin x-x(i)
   		[ b, r ] = deconv( c, [ 1 -x( i ) ] );
   		% calcul p=(x(i)-x(1))...(x(i)-x(i-1))(x(i)-x(i+1))(x(i)-x(n))
   		z = x( i ) - x;
   		z( i ) = 1;
   		p = prod( z );
                a( 1 : n ) = a( 1 : n ) + y( i ) * b( 1 : n )' / p;
	endfor
	a = fliplr(a')';
endfunction

